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A simple model to calculate a-decay Hindrance Factors is presented. Using deformation values 
obtained from PES calculations as the only input, Hindrance Factors for the a-decay of Rn- and 
Po-isotopes are calculated. It is found that the intrinsic structure around the Fermi surface deter- 
mined by the deformed mean field plays an important role in determining the hindrance of a-decay. 
The fair agreement between experimental and theoretical Hindrance Factors suggest that the wave 
function obtained from the energy minima of the PES calculations contains an important part of the 
correlations that play a role for the a-decay. The calculated HE that emerges from these calculations 
render a different interpretation than the commonly assumed n-particle n-hole picture. 

PACS numbers: 23.60.-|-ea decay 27.80.-|-wl90 < A < 219 



I. INTRODUCTION: 

The recent experimental exploration of the neutron de- 
ficient Pb-region has given new insight into the low lying 
excitations of atomic nuclei^]. In particular, the discov- 
ery of the triplet spin 0+ states as the lowest excitations 
in ^^^Pb reveals the coupling of single particle motion 
to the collective shape degree of freedom in a striking 
manner|2|. The quantitative description of these excita- 
tions as well as their isotopic and isobaric dependence is 
a real challenge to nuclear structure models. The balance 
of prolate and oblate shapes e.g. depends crucially on the 
iso-spin dependence of the spin orbit interaction as well 
as the particular shape of the nuclear potential^. 

This paper will focus on the a decay as a tool to under- 
stand the microscopic structure related to the different 
shapes in the Pb-region within the frame work of the 
mean field model. To calculate the absolute values of the 
a decay rate is a problem which is still not fully solved. 
Instead, we will investigate the ratio of the Hindrance 
Factors of the a decay to different states. In this man- 
ner, only parts of the absolute decay width needs to be 
taken into account making the problem more tractable. 
We limit the discussion to a decay of the 0+ ground state 
of even-even nuclei to the 0"*" states in the daughter nu- 
clei; AL = transitions. 

The Hindrance Factor is the ratio between the reduced 
a decay widths, labeled S'^, of the ground state-to-ground 
state decay and the reduced a decay width of the decay 
to an excited state in the daughter nucleus 0. The ex- 
perimental reduced a decay widths are calculated using 
a spherical potential barrier. In principle, using a de- 
formed potential barrier could infiuence the calculated 
(5^-value, but the effect is expected to be very small for 
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the zero angular momentum decays treated in this work. 
The Hindrance Factor is independent of the a-decay en- 
ergy and can therefore be used to compare the intrinsic 
mechanisms that induce the a-decay process. This is an 
important property in the study of shape coexisting nu- 
clei, since the decay between states with similar single 
particle structure is thought to take place more easily 
than that between states with different single particle 
structure. In the present work the single particle struc- 
ture is determined by the deformation parameters of the 
universal parameter Woods-Saxon potential, see ref. |^. 
A study of the Hindrance Factors can therefore help in 
establishing the deformation of the ground state and the 
excited states of the daughter nucleus if the deformation 
of the decaying nucleus is known or vice versa. 



Earlier work on Hindrance Factors have used a two- 
level mixing model based on the spherical shell model 
0. In this model the intrinsic decay rates have been 
related to each other through experimentally measured 
Hindrance Factors |3j 13 ■ ^ microscopic model based on 
RPA calculations has been presented in ref. 0. In the 
present work we start from the mean field wave func- 
tion as obtained from Potential Energy Surface (PES) 
calculations. In this way, the deformation of each state 
is derived microscopically and is not a free parameter of 
the calculations. 



Section 1 of this paper gives an introduction to the 
subject of Hindrance Factors. Section 2 describes the 
model used in our calculations. In section 3 we show 
the dependency of Hindrance Factors on properties such 
as deformation, pairing correlations and particle num- 
ber. Section 4 is a comparison between the nparticle ~ 
nhole {np — nh) spherical shell model picture and the 
mean-field approach applied in our model. In section 5 
we present and discuss results from calculations using the 
model described in section 2. 
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II. MODEL: 

The a-particle decay width can be factorized into two 
parts [Tol |. 



TiR) = PiR)^\m)\', 



(1) 



where P{R) is the Coulomb penetrabihty at the distance 
R between the center of mass and the a-particle, and 
^\T{R)\'^ is the reduced decay width Here T{R) is 
the a cluster formation amplitude. 

The Hindrance Factor is given by the ratio of the abso- 
lute value squared of the a-particle formation amplitudes 
for the decays from the ground state of the mother nu- 
cleus to the ground state (gs) and the excited state (es) 
of the daughter nucleus. This can be written as 



\T,{R;gs^ gs)\^ 
\T;{R-gs^es)\^ 



(2) 



where the ground state to ground state decay is taken 
to be non hindered, i. e. the Hindrance Factor is equal 
to one. Further, the formation amplitude for a decay 
from state B in the mother nucleus to the state A in the 
daughter nucleus can be written as 0| 



If we rewrite the equation |21 in Fock space we obtain 

Ti{R)^ {<P^;-9A\YCm-^B)- (4) 

By using Thoulcss theorem, see , it is possible to con- 
nect two HFB vacuum states with different deformations. 
Therefore the mother nucleus can be expanded in terms 
of two pair excitations of the daughter nucleus, although 
the deformation of the two states may differ. For a more 
thorough investigation of the problem, see ref. El. An 
expression for the mother nucleus as a four particle exci- 
tation of the daughter nucleus is obtained. 



(5) 



where A:, fc is a pair of protons and / , / is a pair of neutrons. 
The factors B{k) and B{1) are related to the pairing den- 
sities 



B{k)^{^B\ctct\^A), 



(6) 



and can readily be evaluated using the Onishi formula 
yielding the expression 



^Blctc+l'fA) - (^b|*a) i-UXU^' V;^),-„ (7) 



where 



and the Ui and Vi are the matrices of the Bogoliubov 
transform for the daughter and mother nuclei. This gives 
a total expression for the formation amplitudes as 

TiiR) = {<S>^;^A\Yi{R)\J2B{k)B{l)c+c^c+cp,^A) 

kl 

= {^A\^A)Y.B{k)B{l) ($„|y,(i?)|c+c+c+c+).(9) 

kl 

The integral can be solved analytically by expressing the 
wave function of the a-particle in terms of Harmonic Os- 
cillator wave functions and using the Moshinsky brackets 
E3 to transform to a center of mass system. However it 
turns out 01 that for the special case of an a-particle 
with angular momentum all terms yY^'^-^f^'^^^''-^^'' of the 
integral 



A[,r- = ($o|>o(i?)|c+c+44) 



(10) 



arc positive and can as a first approximation be replaced 
by an average value Aqq'' . If we now evaluate the expres- 
sion for the decays between two 0+ states we get 



HF = 



T{R- 0+ ^ 0+) _ Ag-(fl) Y.ki B{k)B{l) 
HR- 0+ ^ 0+) Ago"^(i?) Y.k'1' B{k')B{U) 

EkiB{k)B{l) 
J:,,, B{k')B{U) 



(11) 



U ^(U+Ub + V+Vb), 



(8) 



where the 0^ is the ground state and is the excited 0+ 
state. The physical meaning of equation 1111 is that the 
pair transfer amplitude for protons times that for neu- 
trons is the essential quantity that determines the ease 
with which the mother nucleus decays into a correspond- 
ing state of the daughter nucleus. 



III. DEPENDENCY OF THE 
HINDRANCE-FACTOR ON PAIRING AND 
DEFORMATION: 

The deformation affects the Hindrance Factor in two 
ways. First, the degeneracy of levels at spherical shape is 
lifted. This means that two nuclei with different shapes 
have the levels rearranged and different levels become 
occupied at the Fermi surface. Without pairing inter- 
action this would mean that the overlap between two 
wave functions describing a nucleus at different shapes 
becomes zero as soon as we encounter a level crossing. 

Secondly, the spherical components building up each 
single particle level change smoothly with deformation. 
Therefore the hindrance factors are not constant even in 
absence of level crossings. 

The pairing interaction results in a smoothing of the 
differences in level structure for nuclei with different 
shapes. This effect only occurs in the vicinity of the 
Fermi surface, in an energy window that is determined 
by the pairing strength. As a large part of this work fo- 
cuses on the Pb-rcgion we have the further complication 
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FIG. 1: Hindrance Factor as a function of P2 deformation of 
^^°Pb for tfie decay ^^'^Po ^^°Pb. The /Ja deformation of 
^^^Pb is set to 0.1 and both nuclei are assumed to have the 
triaxiality parameter 7 = 0. The result is normalized so that 
the HF of the most favorable decay is set to one. Note that 
the minimum occurs at /32 = 0.12 in the daughter nucleus for 
this specific decay. The increase of hindrance for the spherical 
shape is a consequence of the decrease in pairing correlations 
due to the 82 shell gap. 



HF 




-60 -50 -40 -30 -20 -10 

7 deformation 



FIG. 2: Hindrance Factor as a function of 7 deformation for 
the decay ^^*Po ""Pfe. The 7 deformation of ^^*Po is 
set to 0, both nuclei are set to have /?2 = 0.2, The result is 
normalized so that the HF of the most favorable decay is set 
to one. The inset shows the same plot in logarithmic scale. 

of the large proton energy gap between proton numbers 
82 and 84. The energy gap does not allow for a non trivial 
solution to the BCS-equation. Since the vanishing of the 
BCS solution is non physical and in order to mimic the 
correlations that are indeed present in the ground state 
wave function we have increased the pairing strength by 
15% which guarantees a non trivial solution. The odd 
even mass difference we obtain are well within the ex- 
perimental ones, see also the discussion in ref This 
procedure is necessary in order to compare decays to dif- 
ferent states in a meaningful manner. 



Figure n shows that the HF is only weakly dependent 
on deformation in an interval from P2 ~ —0.1 to (32 ~ 0.2. 
By comparing to Fig. 13 one notices that this is the range 
in deformation, where there are no level crossings for pro- 
ton number Z=82. However, beyond these values in (32 
deformation, we encounter one level crossing at oblate 
shape and several at prolate shape. This is reflected in 
Fig. ^ where the increase in HF is considerably stronger 
at the prolate side compared to the oblate side. Note the 
increase around spherical shape due to reduced pairing 
correlations as a consequence of the Z=82 shell gap. 

In fig. 121 the dependence on the 7 deformation param- 
eter is shown. Similarly to fig. ^ in a rather broad range 
of deformation, from 7 = to 7 = 30 the increase in 
the HF is rather modest. However, for larger difference 
in 7 between mother and daughter nucleus, the HF in- 
creases steeply. In fact, the increase in HF is exponential, 
as shown in the inset. In the calculations of Fig. [21 we 
choose the mother nucleus to be prolate. When starting 
from an oblate deformation for the mother nucleus, we 
obtain a similar curve, but reflected at 7 = —30. Our cal- 
culations imply, that the decay from oblate to prolate is 
essentially forbidden. However, we also investigated the 
case when the mother has maximum triaxial deforma- 
tion, i.e. a 7 value of -30. Then the hindrance increases 
only with a factor of « 4 when the daughter is prolate or 
oblate, respectively. Since true wave function has a cer- 
tain spread in 7 and (32 , the decay from prolate to oblate 
will proceed via the tail of the wave function. This is of 
course beyond the mean field prescription. A soft poten- 
tial in 7, resulting in a large spread of the wave function 
will result in a pronounced decrease of the HF. 



IV. COMPARISON BETWEEN THE 
PARTICLE-HOLE PICTURE AND THE 
MEAN-FIELD PICTURE 

Shape coexistence has been described in the shell 
model in terms of specific particle-hole excitations with 
respect to the spherical core, see In the Pb-region, 
the 2p-2h, 4p-4h or 6p-6h proton excitation across the 
magical shell gap induce deformation that can result in 
oblate or prolate shapes. The spherical ground states 
in the Polonium isotopes are characterized by the corre- 
sponding 2p-0h states and the 0p-2h states correspond 
to the spherical Hg-isotopes. The 2p-2h excitations are 
very costly in energy. On the other hand, the energy 
gain due to pairing interaction and iso-scalar quadrupole- 
quadrupole can partly balance the energy loss and result 
in low lying deformed states. Hence, in a qualitative 
fashion, one has been able to account for the parabolic 
dependence as a function of neutron number of the de- 
formed states in the neutron deficient Pb region as well as 
other regions of the nuclear chart, see e.g. the discussion 
in Ref. 

The shell model description can be linked to the de- 
formed Nilsson (or Hartree-Fock) scheme, by associating 
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the particle hole excitations to the specific level cross- 
ings, present in the Nilsson diagram, which pin down the 
underlying structure of deformed states, see ^3- This 
connection is very valuable, since it enables the compari- 
son of the underlying structure for the two basic nuclear 
structure models. However, two aspects need to be clar- 
ified in this respect: i) within the shell model, there is 
no real distinction between oblate and prolate shapes. In 
the Sn-region e.g., 2p-2h excitations corresponds to pro- 
late shapes in the mean field model, whereas the same 
excitation in the Pb region results in oblate shapes, ii) 
care has to be given to the labeling of the np-nh excita- 
tions. The 2p-2h — (si/2)~^ configuration in the 
Pb-region e.g., does not agree at all with the microscopic 
structure of the corresponding oblate states in the mean- 
field. In fact, the '5^/2' hole state at oblate deformation 
is dominated by the /2 according to the Woods-Saxon 
model, and has only a smaller contribution from the si/2, 
see ref. [T^ . The configuration mixing that is induced by 
deformation is very pronounced for this particular case. 
This is of course essential for a microscopic understand- 
ing of the deformed states, since the quadrupole moment 
of the c?5/2 hole is very different from that of the Si^2- 
Realistic shell model calculations in the Pb region with 
a broken proton and neutron core are far from being fea- 
sible, but structure information from the deformed mean 
field can be used as an important input to shell model 
calculations. 
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FIG. 3: Nilsson diagram for protons in the Z=82 region. Ar- 
eas of interest from a level- crossing perspective are indicated 
both for an oblate and for a prolate deformation. 

Allowed and forbidden decays in the np-nh model can 
be illustrated with figure 2. of ref. The neutrons are 
considered as spectators not influencing the hindrance 
of the decay. Therefore the hindrance of the a-decay is 
discussed in terms of two particle (proton) removal. An 



allowed decay is characterized by the removal of two par- 
ticles from the mother nucleus forming an existing state 
in the daughter nucleus. The mother nucleus is in its 
ground state considered to be a 2-particle state of the 
daughter nucleus. There are two allowed decays from a 
spherical ground state; removal of the two valence parti- 
cles gives a decay to the spherical state in the daughter 
nucleus; removal of two particles from the core forms a 
2p-2h state in the daughter nucleus corresponding to an 
oblate shape. The decay to a prolate 4p-4h state is for- 
bidden since it is a two step process. First two particles 
have to be removed from the core, and after that another 
2p-2h excitation has to take place giving a final 4p-4h 
state. This means that the only allowed decay from a 
prolate state in the mother nucleus is to a prolate state 
in the daughter nucleus. One may also state that the a- 
decay can only remove particles, not holes, i.e. a 4p-4h 
configuration will remain at least a 4h configuration in an 
allowed decay. The different decays are attributed with a 
hindrance taken from the assumed pure (unmixed) decay 
between ^^^Po and Q. Once these pure states are 

mixed, a formula for the Hindrance Factors can easily be 
obtained. The decay properties of these mixed states can 
be used to determine the mixing parameters. 

In the mean-field picture the different 0^ states are 
viewed as 0-quasi particle states with different deforma- 
tions, described by the standard deformation parame- 
ters 7,/32 and /?4, see ref. 0. From figure Q one can 
conclude that the underlying single particle structure is 
changing with deformation, due to the presence of dis- 
tinct level crossings at oblate and prolate shapes. This is 
reflected by the slope of the hindrance factor as a func- 
tion of fi^- However, level crossings are not restricted to 
spherical shapes but a generic feature of a many-body 
fermionic system. In order to obtain a deeper under- 
standing of the HF with respect to the shell model like 
description, we also investigated the neutron and proton 
contributions to the HF. 

20 I , , , , , , 
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FIG. 4: Contribution from neutrons to Hindrance Factor in 
figure 

According to Eqs. I6I1UI the HF are obtained as the 
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FIG. 5: Contribution from protons to Hindrance Factor in 
figure Q 



product of proton and neutron contributions, respec- 
tively. Therefore, one can easily decompose the HF ac- 
cordingly, see Figs. ^ and |S1 Interestingly, the neutron 
contribution to the hindrance factor is of similar strength 
as the proton part. Minor differences are seen at spheri- 
cal shape and at the largest deformations. Since the neu- 
trons in general are treated as spectators, see e.g. Ref. 
our results suggest that this is not the case. Therefore, 
in spite of the appealing clarity of the intruder picture in 
terms of the np-nh excitations with respect to the magic 
shell gap at Z=82, there is little evidence that the Hin- 
drance Factor is dominated by the change of the proton 
configuration. In other words, even if the proton Fermi 
level would be placed in the middle of the shell, as it is 
the case for the neutrons, the calculated HF would be 
similar. One may conclude, that the similarity of proton 
and neutron contribution to the HF indeed indicates that 
this is a generic feature of the nuclear spectrum, and not 
so sensitive to the details of the shell structure. 



V. RESULTS AND CONCLUSIONS 

The deformation values used for the calculations of the 
hindrance factors in eouation llll are determined from the 
minima of the PESs, see figures El [7| |H1 Extensive PES 
calculations of the Po isotopes have been presented ear- 
lier in ref. [3 • We start with the decays of selected Po 
and Rn nuclei, which are calculated to have well sepa- 
rated minima. The results of the calculations arc pre- 
sented in table 1. Given the simplicity of the model, the 
calculations agree surprisingly well with the data. The 
results for the decays of the heavier Po to Pb are better 
reproduced than that of the Rn to Po decays, possibly 
indicating that the deformation of the oblate states in the 
Po isotopes is slightly underestimated in the PES calcu- 
lations. Still, the trend is nicely reproduced. For the 
lightest Po isotope, ^^^Po, the prolate state is calculated 
to have the lowest energy, which is why the hindrance 
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FIG. 6: Potential Energy Surfaces of: a) "^Po, b)"^Po, 
c) Po and d) ^®*Po, with energy minima labeled with dia- 
monds. The energy spacing of the contour lines is 50 keV. The 
X-axis represents /32 cos(7 -I- 30) and the y-axis jii sin(7 + 30) 
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FIG. 7: Potential Energy Surfaces of: a) ^****P6, hf^^Ph, 
c) ^^^P6 and d) ""Pfe, with energy minima marked with dia- 
monds. The energy spacing of the contour lines is 50 keV. The 
X-axis represents (5-2. cos(7 -|- 30) and the y-axis /32 sin(7 + 30) 



of the decay to the excited (prolate) state in ^^^Pb is 
less than one. Also for this case, we reproduce the trend 
of this change, but in absolute numbers the calculated 
Hindrance Factor is three orders of magnitude off. 

A remark is necessary concerning the important role 
played by the small quadrupole deformation, f32 0.1, 
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FIG. 8: Potential Energy Surfaces of; a) Rn, hf°°Rn with 
energy minima marked with diamonds. The energy spacing of 
the contour hnes is 50 keV. The x-axis represents /32 cos(7 + 
30) and the y-axis /32 sin(7 + 30) 



TABLE L Calculated Hindrance Factors compared to exper- 
imental ones. 



Mother Nucleus HFrhe. HFex; 





2.82 


3.2(5)" 


196 


2.38 


2.6(2)" 


194p^ 


2.41 


1.2(2)" 


188 


1.9E-5 


0.08(3)" 
19(6)' 




5.46 


i?n 


10.4 


85(7)^= 



of the Po isotopes, sec fig. El This slight deforma- 
tion, in combination with the reduced pairing correla- 
tions at spherical shape, is essential to reproduce the 
similar strengths of the decay to the spherical and oblate 
states, respectively, in the P&-isotopes. If the Po iso- 
topes would be spherical in the mean field calculations 
the decay from spherical to oblate would be strongly sup- 
pressed. This is demonstrated in Fig. |51 showing that for 
spherical shapes, the decay to either oblate or prolate is 
essentially forbidden. This is also at variance with the 
simple np-nh shell model picture where the decay from 
spherical to oblate is allowed. 

The steep increase in hindrance for the decay of ^^^Rn 
has earlier been attributed to a phase transition, see 
where the structure has changed from the decay of the 
heavier isotope ^'^^Rn. In our calculations the increase 
of hindrance is an effect of the increased deformation of 
the excited state in ^^^Po compared to that of the heavier 
isotope ^^^Po. Since we are in a region of a level crossing, 
see figures n and 13 the calculations are very sensitive to 
small deformation changes as they give a large change in 
single particle structure. This is also in good agreement 
with the PES calculations presented in figures El and (HI 
where only a slight difference in deformation between the 
two Radon isotopes doubles the hindrance of the decay. 
One notices from figuresOandnthat only slight changes 
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l3'2 deformation 

FIG. 9: Hindrance Factor as a function of /32 deformation of 
^^Opfo for the decay ^^-^Po ^^°Pb. The [32 deformation of 
^^*Pfo is set to 0.0 and both nuclei are assumed to have the 
triaxiality parameter 7 = 0. The result is normalized so that 
the HF of the most favorable decay is set to one. 



TABLE II: Calculated Hindrance Factors compared to ex- 
perimental ones. The normalization is such that the decay 
between the spherical states are set to 1 as this proved to be 
the best way of comparing the calculated and experimental 
data. The states are labeled as 1. spherical=s, 2. oblate=o, 
3. prolate=p. 



Mother Nucleus State 


in Dauj 


;hter HFrhe. 


HFexp. " 




0+ 


1.0 


1.0 


''^Po, 


0+ 


1.9 


0.56(7) 


'''^POs 


0+ 


2.8 


>50 


^""^POo 


0+ 


57 


1.0 


''^POo 


0+ 


0.20 


0.56(7) 


'^'POo 


0+ 


8.4 


>50 


^'^OPO, 


0+ 


1.0 


1.0 


^^opo. 


0+ 


1.3 


0.57(12) 


^^opo. 


0+ 


3.5 


2.4(9) 


'^"POo 


0+ 


44 


1.0 


^•'^POo 


0+ 


0.17 


0.57(12) 


^'^''POo 


0+ 


63 


2.4(9) 


'''^POp 


0+ 


300 


1.0 


'"""POp 


0+ 


3.4 


0.57(12) 


'^°POp 


0+ 


0.20 


2.4(9) 



in the equilibrium deformation are necessary in order to 
reproduce experimental data. 

We also calculated the Hindrance Factors for the the 
decays ^'^^Po -> ^^^Pb and ^^"Po ^ ^^^Pb. For these 
decays the different minima in the PES of both mother 
and daughter nuclei are less well separated. In general, 
one expects the mean field to give a proper description 
only when the energy minima are sufficiently deep, and 
the spreading of the wave function can be neglected. Cer- 
tainly this is not the case in the light Po and Pb isotopes. 
However, for pedagogical reasons and also to indicate the 
limitations of the present approach, we compare the cal- 
culated and experimental values in table 2. 

An interesting feature in table 2. is the coupling be- 
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tween states of difFercnt deformations which is not in 
agreement with the particle hole picture, see ref. @. 
Note also the influence of triaxiality visible from the dras- 
tic change in hindrance for the decays oblate to prolate in 
the two different Po isotopes. The triaxiality parameter 
7 is -13.9° for ^^^Pb and -3.7 for ^^^Pb. The separa- 
tion of the minima in deformation and energy obtained 
from the PES calculations is now the limiting factor of 
calculations. Creating a proper collective wave function 
by means of the Generator Coordinate Method (GCM) 
e.g. is a future project that may enable us to take into 
account the mixing of the wave functions at different de- 
formation, resulting in meaningful calculations of all a- 
decaying shape coexisting nuclei. GCM calculations has 
been performed for the mass region, see ref. , but has 
been restricted to axially deformed states. 

To conclude, a simple model to calculate the Hindrance 
Factors for alpha decay is presented. Both the energy 
minima and the hindrance factor calculations are per- 
formed using the same Woods-Saxon potential with uni- 



versal parameters. Deformation values are taken from 
microscopic calculations and not free parameters. In this 
respect, one may view the alpha decay as an important 
probe to the mean field wave function. Our calculations 
show that the HF factors depend only modestly on defor- 
mation changes, as long as these are not large. However, 
as soon as the underlying s.p. structure is changed due 
to a level crossing e.g., there appears an exponential in- 
crease in the HF with deformation. Our calculations re- 
produce the experimental trends for the cases where one 
deals with clearly separated minima in the PESs. This 
indicates that the mean field wave function contain an 
essential part of the correlations that are responsible for 
the formation of an alpha particle. In this context it is 
important to point out the role of the calculated small 
deformation values in the Po- and Rn-region which arc 
normally treated as spherical. Our calculations also re- 
veal that both protons and neutrons contribute to the 
HF in a similar fashion, which is in variance with the 
simple shell model like description. 
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